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Abstract 


Using the Dirac-Frenkel time-dependent variational principle and the multi- 
Dz Ansatz we have simulated the process of singlet fission. The singlet fission 
process has seen a resurgence in the recent years because of its ability to by- 
pass the Shockley-Queisser limit of photovoltaics. We have studied the effect of 
phonon frequencies associated with diagonal and off-diagonal coupling (diagonal 
and off-diagonal coupling frequencies) and the role of the Huang-Rhys factor. In 
the case of high off-diagonal coupling frequencies, the efficiency is independent 
of the Huang-Rhys factor, while in the case of the low off-diagonal coupling fre- 
quencies, we obtain a high-efficiency channel. We have further studied the time 
scales involved in the process for different diagonal and off-diagonal coupling 
frequencies. We intend to provide a frequency-efficiency analysis to help future 
development of singlet fission materials. 
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1 Introduction 


Solar energy is readily available, clean and a safe source of energy. While there is 
plenty of energy, harvesting it efficiently has been a challenge for a long time. The 
Shockley-Queisser limit, which gives a maximum efficiency, puts a constraint on all 
photovoltaics [1]. One of the ways to go around the limit is to increase the quantum 
yield by generating more than one exciton per photon. This has been demonstrated 
in materials capable of singlet fission (SF) by Hanna et. al [2]. 

SF is a spin-allowed conversion process in molecules and molecular aggregates, 
in which a spin-singlet exciton generated by irradiation splits into two spin-triplet 
excitons [1]. Singh et al., in 1965, first used the SF process to explain the unusual 
photophysics in anthracene crystals [3]. In this process, So, the electronic ground 
state and S 1, the lowest singlet excited state follow the kinetic model Sg + S$; > 
TT — T+T, where T is the molecular triplet state and TT is a doubly excited pair 
of spin-correlated triplets which has an overall singlet spin. The TT state is called 
a correlated triplet pair state and is considered as a dark state because it cannot be 
optically populated from the ground state. 


SF has received a great deal of attention again since an external quantum ef- 
ficiency above 100 % has been realized in a SF-based organic photovoltaic cell by 
Congreve et. al [4]. In order to devise more rational design of photovoltaic sys- 
tems, many experimental [5, 6, 7, 8, 9] and theoretical [10, 11, 12, 13] efforts have 
been devoted to gain further understanding of the SF process at the molecular level. 
However, the current limited understanding of detailed SF mechanisms hinders the 
design of versatile SF materials. In particular, a unified treatment of phonon effects 
remains elusive. 

The SF process has been found to be influenced by various factors, in which 
molecular vibrations were predicted to significantly affect the SF efficiency because 
the exciton-phonon coupling served as the origin of the exciton relaxation [10, 11]. 
Recently, phonon modes coupled to electronic excitations have been shown to play 
a crucial role in the SF process according to ultrafast spectroscopic measurements 
in SF materials [6, 7, 14]. For example, efficient fission in pentacene derivatives 
(15, 16, 17] and crystalline tetracene [18, 19] has been facilitated by high-frequency 
phonon modes because of resonances between vibrational modes and energy split- 
tings of electronic states. It is thus essential to explicitly explore the impacts of the 
vibrationally induced fluctuations on the fission dynamics through quantum dynam- 
ics calculations. 

There exist two types of exciton-phonon coupling, i.e., diagonal coupling and 
off-diagonal coupling. The diagonal coupling term represents fluctuations of the 
energy gap between the optically-allowed state S; and forbidden state TT induced by 
intramolecular vibrations. The off-diagonal coupling term describes the fluctuations 
in excitonic coupling induced by intramolecular and intermolecular vibrations. In 
covalent tetracene dimers, using quantum mechanical calculations, was uncovered 
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that high-frequency intramolecular vibrations give rise to strong diagonal and off- 
diagonal coupling [20]. The two types of coupling have also been found to be tunable 
in novel SF materials by changing linkers and by engineering the dihedral angle 
between the chromophore units and the linker [21, 22, 23]. However, there is a lack 
of discussion in the literature on detailed SF mechanisms under the influence of 
simultaneous diagonal and off-diagonal exciton-phonon coupling. 

Impacts of diagonal and off-diagonal coupling on exciton dynamics in organic 
crystals have been investigated by Zhao and coworkers, using a refined trial state, 
the multiple Davydov D2 Ansatz [24, 25], to accurately treat dynamics of the Holstein 
model [26, 27]. Recently, influences of diagonal and off-diagonal coupling have also 
been probed in the intramolecular SF model using our variational approach [15, 28]. 
We have found that both diagonal and off-diagonal coupling can aid efficient SF 
if excitonic coupling is weak. But it remains an open question on whether those 
conclusions are valid if the Huang-Rhys factor changes. 

This work is a continuation of the previous investigation, answering the question 
by using a different set of Huang-Rhys factors. The method used in this work is still 
within the framework of the Dirac-Frenkel time-dependent variation. Here, we will 
only briefly explain the adopted Hamiltonian: 


HA = Fie a Hath + Haek-yerh (1.1) 


The first term of H is an electronically diabatic Hamiltonian for |g), |S;) and |TT) 


Hoys= > engln)(nl + SY) Yo Fini) (n, (1.2) 


n=S51,TT m=S1,TT nxm 


where €ng is the Franck-Condon energy associated with the electronic transition from 
|g) to |n), and Js, rr is the strength of the interstate coupling between $; and TT. 
Js,,rr includes the contribution of the direct coupling between S; and TT based on 
the two electron integrals {1, 10] as well as that of the effective coupling created by 
quantum mixing of the charge transfer state and the electronic states. The second 
term of H represents the bath Hamiltonian Hyatn, and is given by 


bath = a hugdl bg, (1.3) 
q 
where w, indicates the frequency of the qg-th mode of the bath with the creation 
(annihilation) operator, bf (b,). The final term represents the system-bath coupling, 
and is given by 


Asys—path = Len=81,TT [72) (r|Rurgaq (bi, + bq) 
oP See ie ore |) (| fiwrgeq(b + bg); (1.4) 


where gq and cy are the diagonal and off-diagonal exciton-phonon coupling strengths 
between the system and qg-th mode, respectively. The spectral density Ja(w) (a = 
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diag,o.d.) is a useful measure for characterizing various forms of exciton-phonon 
couplings, and can be evaluated in terms of gg (cq) as 


dnl) = SO hu k70(w — Wg), (1.5) 


where kg = gq for diagnonal coupling and cg for off-diagonal coupling. In this 
study, we model spectral densities using underdamped Brownian oscillators with 
the diagonal coupling Huang-Rhys factor, S;, = Amg/(hwaiag), and an off-diagonal 
coupling Huang-Rhys factor, $04 = \°-4-/(fhw,.g,), such that 


Ay ate 
Jo(w) = Woah ea? (1.6) 
where We, is the vibrational coupling frequency and Yq is the vibrational relaxation 
rate corresponding to diagonal and off-diagonal coupling for the respective @ = 
diag,o.d. Amn represents the reorganization energy associated with the transition 
from |m) to |n), and A°:*- is the amplitude of fluctuations in the interstate coupling 
between |m) and |n). Detailed descriptions for the derivation of the equations of 


motion for the variational parameters are available in Ref. [28]. 


2 Results and Discussion 


Using the multi-D2 Ansatz we have simulated the system to a high degree of accuracy. 
In this work we will present numerical results for the SF process. A new set of Huang- 
Rhys factors, other than that in Ref. [28], is used to study the effect of phonon 
mode frequencies associated with diagonal and off-diagonal coupling (diagonal and 
off-diagonal coupling frequencies, hwajag and hwo.a,). The Huang-Rhys factor for 5; 
is Ss, = 0.7, estimated from fitting measured absorption spectra of acene derivatives 
with a spectroscopic model [29]. Since the reorganization energy of TT is several 
times larger than that of S; in pentacene derivatives and tetracenes [12], the Huang- 
Rhys factor for TT is set to be Spr = 3Sg, instead of Spr = 2Sg, in Ref. [15, 28]. 
The Huang-Rhys factor for off-diagonal coupling Soci is chosen to be 0.2 instead 
of 0.1 [28] to examine the SF process. 

To understand the effect of different parameters and identify the energies and 
phonon frequencies on the process, we have graphically represented different scenar- 
ios. First we consider the relatively simpler cases where the excitonic coupling is 
weak (Js, rr = 20 meV) and the energy gap is large (€s, . — err,g = 100 meV) while 
staying within the range for the values observed for pentacene derivatives. In the 
absence of the exciton-phonon interaction, the oscillations are completely localised, 
and the amplitude locked at 0.13. When the vibrations are present, we have a differ- 
ent, more interesting process taking place. As we can see in case for both diagonal 
and off-diagonal coupling in Figures la and 1c, there is a very rapid drop in the 
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Figure 1: Time evolution of the singlet population for the case of weak excitonic 
coupling Js, rr = 20 meV, €s,,¢ — ert,g = 100 meV in the first column and for the 
case of strong excitonic coupling Js, rr = 80 meV, €s,.¢ — ert.g = 30 meV in the 
second column. In each subplot, the green line refers to the case without coupling. 
The blue lines correspond to the diagonal coupling frequency of hwgiag = 65 meV in 
the upper panels and the off-diagonal coupling frequency of hwo.q¢, = 75 meV in the 
lower panels. 
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Figure 2: Snapshots of singlet population as functions of hwgiag and hwo, at the 
time of (a) 0 fs, (b) 250 fs, (c) 500 fs and (d) 2000 fs. The other parameters are 
Js,,7TT = 20 meV, €s, ¢ — ert.g = 100 meV, and Vaiae = V4 =1ps. 


—hwiiag = 25meV —Tiudging = 25meV 
—hw4iag = 45meV —fwaiag = 45meV 
0.8 —hwaiag = 25meV —hwaing = 65meV | —hwaiag = 65meV 
oO —hwaiag = 45meV —hwaiag = 85meV —hiwhiag = 85MeV 
© 0.6 —hwiag = 65meV 
s hwiiag = 85meV 
Ss 
2 04 
fo} 
a 
0.2 
(a) hwoa, = 20meV (b) hwo.a, = 60meV 
1 eae eee 
ae . Nal Wamp gpadalar tanreenrirengea aE | 
0.8 
ao —hwiag = 105meV 
3 06 —hwyiag = 125meV 
5 5 —hwiag = 145meV —fwating = 105meV 
= 0.4 —hredtiog = 165MEV | Auyigg = 125meV Rwiag = 165meV 
a : —hwaiag = 145meV 
7 —hwuiag = 165meV 
0.2 : 
(d) fwd. = 20meV (ec) hwoa. =60meV ‘| (f) Awy.a, = 100meV 


% 500 1000 1500 0 500 1000 1500 0 500 1000 1500 2000 
t(fs) t(fs) 


Figure 3: Time evolution of singlet population. First, second, and third columns 
correspond to phonon modes of hwo.q, = 20 meV, 60 meV, and 100 meV that are 
off-diagonally coupled to the S; and TT states, respectively. The phonon modes 
diagonally coupled are hwdiag = 25,45,65 and 85 meV (upper panels) and hwgiag = 
105, 125,145 and 165 meV (lower pannels). The other parameters are same as Fig. 2. 
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population of S; which is converted to the TT state. Here we have shown that for 
hwaiag = 65 meV in Figure la and and hwo.q, = 75 meV in Figure 1c, there is very 
efficient SF happening. 

On the other hand, the strong excitonic coupling leads to different behaviour. 
Due to the strong excitonic coupling and the absence of the phonons, the population 
of 5; shows a pure Rabi oscillation. This is because the stronger coupling facilitates 
the oscillation between S; and TT. Thus, we have a state where no effective SF takes 
place even at longer durations. When the exciton-phonon coupling is introduced to 
this system, we find that for hwgiag = 65 meV we again get an efficient SF happening 
as the population is below 0.5. But with only off-diagonal coupling, the process is 
inefficient. 

Due to the fact the process is not always efficient, and it depends on both types 
of phonon frequencies, diagonal and off-diagonal, coupled to the excitonic states, we 
need to study the effect of both types on a single system. To do this we considered 
the case of Js, rr = 20 meV and €g, 5 — err.,g = 100 meV as where the population 
is plotted on a contour map with the diagonal and off-diagonal coupling frequencies 
as the x and y axis respectively (Figure 2). As we can see from the figure, diagonal 
coupling frequencies give efficient SF only in a narrow band for low coupling strength 
which is centered at hwgiag = 65 meV. Comparing with our previous results, we 
notice that the narrow channel is always present and only a shift in the location along 
the hwqiag scale happens. The location is affected by the Huang-Rhys factor but the 
overall structure remains consistent. This implies that, depending on the Huang- 
Rhys factor, there will be a certain frequency of the diagonally coupled phonons 
which will give efficient SF. The combination of high off-diagonal and low diagonal 
coupling frequencies will always give efficient SF. 

To complete the analysis and understand the time-scales involved in the process, 
we have calculated the effect of diagonal coupling frequencies in three different ranges 
for off-diagonal coupling frequencies (Figure 3). For this purpose we have selected 
three different regions (weak, medium and strong) and use three representative val- 
ues of hiWo.q.= 20, 60 and 100 meV. Also, we have divided the cases of diagonal 
coupling frequencies into two different sections corresponding to weak and strong 
couping. In case of the weak off-diagonal coupling frequencies we observe that with 
the exception of hwgjag = 65 meV, which also depends on the Huang-Rhys factor 
chosen, the process is always inefficient. This is in agreement with the calculations 
made earlier [28]. 

In case of the medium off-diagonal coupling frequencies, i.e. hwo.g, = 60 meV, 
we find that the case of hwgiag = 65 meV is efficient even at very short duration, 
while the cases with other coupling strengths require some time to attain efficiency. 
Increasing coupling strength from hwgjag = 65 meV results in a small increase in the 
efficiency followed by a rapid decrease in efficiency. 

For high-frequency off-diagonal coupling phonon mode and low-frequency diago- 
nal coupling phonon mode we find the SF process to be very efficient and the system 
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reaches the TT state in a very short time. It is also important to notice that in the 
long-time range this is the most efficient region available as the singlet population 
drops to approximately 0.1. The efficient fission for the case of hwo.g. = 100 meV is 
different from those for hwo,q, = 20 and 60 meV and gives a unilateral decrease in 
the efficiency with the increase in the hwgjag. 


3 Conclusions 


In this work, we have attempted to characterise the SF process on the basis of its 
dependence on the phonon vibrations corresponding to the diagonal and off-diagonal 
coupling (which are given by w¢jag and w,4,) and the excitonic coupling strength 
between the 5; and TT states (which is given by Js, rr). Using the Dirac-Frenkel 
time-dependent variational principle and the multi-D2z Ansatz, a superposition of the 
usual Davydov Dz states, we have simulated the system preserving accuracy. We 
have presented evidence that both diagonal and off-diagonal coupling frequencies 
play an important role in the SF process. It was also shown that the Huang-Rhys fac- 
tor can significantly affect the interplay between the exciton and phonons especially 
in the case of weak off-diagonal coupling as the variation of the Huang-Rhys factor 
can alter the location of the efficient fission channel. In case of the high-frequency 
off-diagonal coupling modes the efficiency is independent of the Huang-Rhys factor 
and decreases with an increase in the diagonal coupling frequencies. Our work here 
provides a guideline to facilitate the development of novel SF materials. 
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